Thermodynamic properties of Holstein polarons and the effects of disorder 
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The ground state and finite temperature properties of polarons are studied considering a two-site 
and a four-site Holstein model by exact diagonalization of the Hamiltonian. The kinetic energy, 
Drude weight, correlation functions involving charge and lattice deformations, and the specific heat 
have been evaluated as a function of electron-phonon (e-ph) coupling strength and temperature. The 
effects of site diagonal disorder on the above properties have been investigated. The disorder is found 
to suppress the kinetic energy and the Drude weight, reduces the spatial extension of the polaron, 
and makes the large-to-small polaron crossover smoother. Increasing temperature also plays similar 
role. For strong coupling the kinetic energy arises mainly from the incoherent hopping processes 
owing to the motion of electrons within the polaron and is almost independent of the disorder 
strength. From the coherent and incoherent contributions to the kinetic energy, the temperature 
above which the incoherent part dominates is determined as a function of e-ph coupling strength. 



I. Introduction 



Study of different properties of polarons has been of great 
importance since the evidence of polaronic charge carri- 
ers in many materials of recent interest, viz. high-Tc 
cuprates"'^, CMR-manganites^, biological materials like 
DNA^, etc. which have large technological potential and 
importance. In the simplest Holstein model an electron 
in a narrow tight-binding band interacts locally with dis- 
persionless optical phonons. For large e-ph coupling the 
resultant polaron is a small polaron with high effective 
mass, while for weak coupling it becomes a large polaron 
having a much lower effective mass for a finite adiabatic 
parameter. The crossover from a large to a small polaron 
and the corresponding changes in the polaronic proper- 
ties in the ground state have been studied for the Hol- 
stein model by different groups^"^* using various meth- 
ods to enrich our understanding in this field. However, 
finite temperature study of the properties of polarons and 
the effect of disorder on polaronic properties are few and 
needs more attention. Previously we have studied the 
effect of disorder on some of the polaronic ground state 
properties (i) for a two-site system following a pertur- 
bation method based on a modified Lang-Firsov (MLF) 
phonon basis"'^^ and (ii) for a many-site system following 
a zero-phonon averaging of the MLF-transformed hamil- 
tonian and a real space renormalization group method 
to deal with the disorder^^. The above studies have the 
limitations that they are only for the ground state. The 
first study, though quite accurate, has been carried out 
only for a two-site system, while the latter study gives 
approximate results. In this paper we will consider a 
two-site and a four-site Holstein model and follow an ex- 
act diagonalization method to study the ground state as 
well as finite temperature properties of the polarons and 
the effect of disorder on them. We will mainly study the 
kinetic energy, correlation functions involving charge and 
lattice deformations, Drude wight and the specific heat 
of the systems as a function of e-ph coupling for different 



temperatures and disorder strength. 

The paper is organized as follows. In Section II we have 
developed the formalism for the aforementioned study 
considering the Holstein model. We have presented the 
results and discussions for the two-site Holstein model 
in Section III-A and those for the four-site system in 
Section III-B. The conclusion is given in Section IV. In 
Appendix-A we have shown analytically considering an 
infinite size system and following strong-coupling second 
order perturbation theory that the disorder has weak or 
negligible effect on the kinetic energy for strong coupling. 

II. Formalism 



The Holstein Hamiltonian with site diagonal disorder in 
1-d is given by 

= ^ eicfci - ^(tejci+i + h.c) + guj'^ni{bl + bi) 

i i i 

+ ojY,b\h (1) 

i 

c\ and Ci are the electron creation and annihilation oper- 
ators at the site i, rii (= c\ci) is the number operator, h\ 
and hi are the creation and annihilation operators for the 
phonons corresponding to interatomic vibrations at site i 
and Lo is the phonon frequency. Electronic hopping takes 
place only between the nearest-neighbor sites with hop- 
ping strength t and g denotes the local e-ph coupling. 
The electronic site energy is independent of the site 
i for the ordered case. To study the effect of the site- 
diagonal disorder we would put a different site potential 
at one of the sites of the 2- or 4-site system. Spin index 
is not used for the electron, because a single polaron case 
has been studied here. 

The third and fourth terms of Eq.l represent the 
electron-phonon interaction and phonon harmonic en- 
ergy, respectively. These terms may be written in the 
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momentum space defined by the phonon creation oper- 
ators: 6^= {l/^/N) Y^ible^^-^^ and the corresponding 
annihilation operators, where N is the number of sites in 
the system'"' It can be easily shown that the in-phase 
(q — 0) phonon mode does not couple with the electron 
dynamics but with the total number of electrons of the 
system. The harmonic term of this phonon mode along 
with its interaction with the electron may be separated 
out and written in a diagonal form"'''^. The rest of the 
Hamiltonian involving {N— 1) phonon modes and N elec- 
tronic states (for a single electron problem) are consid- 
ered to construct the cigen basis and matrix elements for 
diagonalization of the matrix. If one considers Up num- 
ber of phonon states per mode then the total number of 
basis states will be utol = Nrip^ . Elimination of the 
in-phase mode, thus, reduces the states of the Hilbert 
space by a factor of Up, which is an advantage for any 
diagonalization procedure. 

For the electron states we use the site space basis, 
which is convenient to take into account of the site disor- 
der. For the phonon states we use the momentum space 
basis so that the in-phase (q = 0) mode may be separated 
out. The Hamiltonian is then diagonalized to obtain the 
eigenstates and the eigenenergies. Thermodynamic ex- 
pectation value of any observable characterized by the 
operator O is then found out by 



-, "Tot 

(0> = ^l](n|0|n)e-'5^" 

"Tot 



-0E„ 



(2) 
(3) 



n=l 



where _E„ is the cigcn energy of the n-th eigenstate 
/3 = l/ksT and T denotes the temperature. In this pa- 
per we are interested in evaluating the kinetic energy, 
static correlation function involving charge and lattice 
deformation, specific heat and the Drude weight. The 
operator corresponding to the kinetic energy is 



(4) 



For the correlation functions involving charge and lattice 
deformations we calculate 



Xm{i) = {ni{bl+^ + bi+m))/2g 



(5) 



which represents the lattice deformation produced at the 
site i+m when the electron is at site i. The Drude weight 
(D„) in units of ire^ for an eigenstate |n) of the polaronic 
system is obtained by introducing a phase factor to the 
hopping matrix element {t —* te^"^) in order to break the 
time reversal symmetry and then finding out the response 
of the break down of the time reversal symmetry to the 
electric current as^'^ 



Dn 



(6) 



where En (</>) is the eigen energy of the n-th eigenstate in 
presence of non zero (p. The thermodynamic expectation 
value of the Drude weight {D) is found out by taking the 
thermal average of Z)„ over all the eigen states 
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(7) 



The specific heat may be expressed in terms of the en- 
ergy fluctuation of the system at a finite temperature T 
as 
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{kBTf 



[{E^) - {Ef 



(8) 



where {E) and {E"^) are the thermal average of the energy 
and the square of the energy respectively. 

III. Results and discussions 



A. Tviro-site system : 

For the 2-site system the electron dynamics is coupled 
only to the out-of-phase (q = tt) phonon mode. The 
Hamiltonian, which has to be considered for numerical 
diagonalization, is^^ 

Hd = ^CiTii - t{c\c2 + c\ci) 

i 

+ wg+{m-n2){d + S)+wSd (9) 

where g+ = g/V^ and d= {b\ — 62) / \/2- We consider the 
basis states cl\{))i,\nfi)ph, where v'=l,2 and nd= 0, 1, 2, 
njp, is the number of phonons in the d-oscillator. We refer 
to this as the bare basis. The Lang-Firsov transformed 
or MLF-transformed basis for the d-oscillators may also 
be used for diagonalization. But the matrix elements of 
the Hamiltonian operator are much simpler in the bare 
basis than those in other (LF or MLF) basis. We find 
that for exact diagonalization study the convergence of 
the results is achieved with much less number of phonon 
states within the bare basis compared to the LF or MLF 
basis, although the MLF basis is best and much better 
than the bare basis for perturbation calculation'^^. 

To show the convergence of the correlation functions 
in the bare basis we have evaluated (niW2), where U2 = 
(&2 + ^2)1 using different values of Up, and plotted it as a 
function of g+ for different values of t and (= €2 — ei) 
in Fig.l. We will refer to Cd for the two-site system as 
the disorder strength since it partly mimics the role of 
disorder in larger systems. It is found that in the range 
< g+ < 3.2 the results obtained for rip = 25 is quite 
accurate. The results for Up = 25, 30 and 40 are indis- 
tinguishable. We also find that for a particular value of 
Up the highest value of g+, up to which the convergence 
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is achieved, almost does not depend on or t (Fig.l) in 
the range of parameters we have studied. 

In Fig. 2 we have plotted the kinetic energy evaluated 
for the parameters t =1 and 6^=1, using different values 
of Up. It is found that the curves are indistinguishable for 
rip > 15 in the range < 3+ < 3.2. It shows that a much 
smaller value of Up is sufficient to achieve the convergence 
in the kinetic energy than that for the convergence in the 
correlation function (m (&2 + &2 )) • the following we will 
present the results for the two-site system for rip =40 or 
50. The latter is used for high temperature and strong 
coupling. 

In Fig. 3 we have presented the variation of the cor- 
relation function Xm{i) = (niUi^m) /^g for i=l,2 and 
m=0,l for the parameters t = 2.1 and = 1. In the 
absence of disorder Xm{i) is independent of i. An in- 
troduction of the disorder (e^ 7^ 0) breaks the transla- 
tional symmetry of the system and xo(l) and Xo(2) be- 
come different such that Xo(l) — Xo(2) increases with the 
increase of the strength of the disorder potential (e^). 
However, it is observed that (niU2) = («2Wi) even in 
the presence of disorder. For a periodic Holstein model 
the on-site correlation involving charge and lattice de- 
formation (riiUi) is always larger than the intersite cor- 
relation (rijWjii). With increasing e-ph coupling (riiUi) 
increases while {niUi±i) decreases indicating a crossover 
from a large polaron to a small polaron. In Fig. 3 it is 
seen that in presence of disorder the correlation (722^2) 
decreases with increasing e-ph coupling and in the strong 
coupling region (n2U2) < {n2Ui). These are very dif- 
ferent from the normal behavior of polarons in ordered 
systems. However, (niUi) and {niu^) show the usual 
characteristics of the polaron crossover with increasing 
e-ph coupling, The above-mentioned observations create 
confusion in distinguishing the large and small polaron 
from the structure of the correlation function Xm{i) in 
the disordered polaronic system. In this situation we de- 
fine average correlation functions 



avg _ T.,Xa{i) 

AO 



avg _ Ei Xlji) 



_EM^)-Xl{i)) 

Ad 



(10) 

(11) 

(12) 



to describe the polaronic behavior for the disordered po- 
laronic system. 

In Fig. 4 we have shown the variations of the kinetic 
energy and X^"^ with g+ for f = 2.1 for different values 
of Cd- It is seen that the kinetic energy is suppressed 
with increasing e-ph coiipling as well as with increasing 
disorder strength. In the intermediate range of coupling 
the kinetic energy shows an exponential suppression. For 
strong coupling the kinetic energy shows a behav- 
ior and is almost independent of the disorder strength, 
as noted previously^^. In the Appendix- A, considering 



an infinite lattice Holstein model we have shown analyt- 
icalUy that for strong coupling the polaronic hopping is 
so small that the kinetic energy arises mainly from in- 
coherent hopping owing to the undirected motion of the 
electron keeping the centre of the polaron fixed. This 
incoherent hopping contribution for large g is almost in- 
dependent of the disorder strength and inversely propor- 
tional to g'^. Fig. 4 shows that with increasing disorder 
strength the correlation x™^ increases, which indicates 
that the size of the polaron becomes smaller, and the 
crossover (from a large to a small polaron) occurs at a 
lower value of g+. For strong coupling, disorder has al- 
most no effect on x™^ as similar to that observed for the 
kinetic energy. It may be noted that increasing disorder 
strength makes the polaron crossover more smoother. 

To examine the effect of temperature on the proper- 
ties of polarons we have evaluated the kinetic energy and 
the correlation function for different temperatures and 
disorder strength. In Figs. 5 and 6 we have plotted 
Xd"^ and the kinetic energy for t = 2.1 as a function 
of for different temperatures. For weak and inter- 
mediate coupling Xd"^ increases, implying that the size 
of the polaron becomes smaller, while the kinetic energy 
is siipprcssed with increasing temperature. For stronger 
coupling {g+ > 1.7) the temperature has an opposite ef- 
fect, i.e., the kinetic energy increases with temperature. 
However, the effect of temperature is small for strong cou- 
pling. At high temperature the variation of the kinetic 
energy with e-ph coupling is much weaker compared to 
that at low temperatures. It is also found that disorder 
has very little effect on both the kinetic energy and the 
correlation function at high temperatures (not shown in 
the figure). 

We have studied the specific heat (Cy) of the ordered as 
well as disordered polaronic system for different temper- 
atures and hopping. In Fig. 7. we have presented the Cy 
of the ordered two-site system as a function of e-ph cou- 
pling for different temperatures for the hopping param- 
eter t=2.1. Similar plots for different disorder strengths 
(ed = 0,0.5 and 1.0), and for t=l and 2.1 are shown in 
Fig.8. 

Fig. 7 shows that in the low temperature regime the 
specific heat shows a peak at intermediate coupling. 
With increasing temperature the peak shifts towards a 
lower value of g+ and then disappears while a dip is 
developed in the intermediate coupling region. At low 
temperatures the specific heat is mainly governed by the 
separation of the ground state and the first excited state 
of the two-site Holstein model. 

It may be noted that the specific heat for a system hav- 
ing only two energy levels with energy separation {AE), 
shows a peak at AE/ksT = 2.58, and Cy is very small 
when AE is far from 2.58A:bT. For the two-site Holstein 
model the energy separation {AE) between two lowest 
eigen energy levels decreases monotonically with 5+ and 
becomes negligibly small at strong coupling-'^* . At very 
low temperatures Cy is very small for weak e-ph coupling, 
as AE/ksT is sufficiently large. Cy attains a maxima 
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at an intermediate coupling, when llE is ^ i.bSksT, 
and becomes very small in the strong coupling region as 
A.E {). At a higher temperature a higher value of AE 
is required to achieve the maximum value in Cy. Since 
AE increases with decreasing for the two-site Holstein 
model, the peak in is obtained at a lower value of g.^ 
for a higher temperature. With increasing temperature 
the higher energy states, in addition to the ground and 
first excited states, have significant contributions to the 
specific heat. This leads to the absence of the peak and 
formation of a dip in Cy. In Fig. 8 a comparison of the 
variation of specific heat as a function of g+ for the hop- 
ping parameters t = 1 and 2.1 is given. In absence of any 
disorder the energy separation A.E increases with the in- 
crease of the hopping parameter t, hence a shift in the 
position of the peak in C„ is obtained at a higher value 
of g+. 

The effect of disorder on the specific heat of a polaron is 
also shown in Fig. 8. In presence of disorder the specific 

heat is suppressed for weak and intermediate coupling. 
However, for strong coupling the C„ is larger compared 
to that for the ordered case. 

B. Four-site system : 

For the four site Holstein model, out of the four phonon 
modes there are three modes with q 7^ which couple 
to the electron dynamics. The four-site Holstein model 
Hamiltonian may be divided into two parts: one part 
containing the electronic terms, harmonic terms of the 
three phonon modes and the interaction of these phonon 
modes with the electron number operators, while the sec- 
ond part gives a diagonal form for the shifted in-phase 
(q = 0) oscillator (see Eqs. 7 and 8 of Ref.6). The 
first part, which contains the non trivial physics of the 
system and cannot be treated exactly by any analytical 
method, is diagonalized numerically. We have done the 
numerical diagonalization with Up—Q per phonon mode 
and checked that this gives fairly accurate results in the 
range of our study presented here. For the disordered 
case we break the translational symmetry by introducing 
a different site-potential at one of the lattice sites while 
keeping the equal site-potential for the rest. 

For the Holstein model the correlation functions satisfy 
a sum rule 

^Xm{i) = ^{niUi+rn)l'^g 
m m 

= {ni) (13) 

In absence of disorder {uiUj) = {rijUi) because of trans- 
lational symmetry and a relation X]i("«'^j)/2.9 = (^j) is 
satisfied. We find that even for the disordered Holstein 
model (riiUj) = (rijUi), though the sites i and j have 
different site potentials, and the relation '^^{niUj) /2g = 
{rij) is also satisfied for the disordered system in addition 
to the sum rule (13). As mentioned previously average 



correlation functions should be used for the disordered 
case to characterize the nature (large or small) of the po- 
laron. Averaging also justifies that the disordered impu- 
rity may occupy any site at random. The average correla- 
tion functions involving charge and lattice deformations 
are: 

avg _ J2iXm{i) 

where m = 0, 1, 2. 

In Fig. 9 we have plotted the kinetic energy (scaled to 
2t) and Xo'^ against g for the ordered case and for disor- 
dered cases with site potentials (1,0,0,0) and (-1,0,0,0), 
where the values of the site potentials (e,) for z = 1, 2, 3, 4 
are shown within the parentheses (all the energies are ex- 
pressed in units of w = 1). It is seen that the polaron 
becomes more localized (in size) with higher value of Xo"^ 
and lower kinetic energy for disordered cases compared 
to the ordered case. For the (-1,0,0,0) case the value 
of Xo"^ is much higher and the kinetic energy is much 
lower than those for the (1,0,0,0) case, because in the 
former case the electron will tend to be trapped at the 
site of the negative potential. This would suppress the 
kinetic energy and favor small polaron formation in pres- 
ence of e-ph coupling. In the same figure we have also 
shown the variations of xi and X2 for the ordered case. 
Xi and X2 represent the lattice deformations produced at 
the nearest- and next-nearest neighbor sites of an elec- 
tron. For small g the values of xi and X2 are appreciable 
indicating that the polaron has spread over the lattice. 
With increasing g, xi and X2 reduce and become very 
small for strong coupling, the corresponding polaron is 
a small one. For the disordered cases the xV^ and xV^ 
behave in the same way (not shown in the figure), but 
their values arc lower than those for the ordered case and 
become insignificant at a lower value of g compared to 
the ordered case. In the strong coupling limit the depen- 
dence of the kinetic energy and the correlation function 
on the disorder strength is very weak. The reasons have 
been explained in Appendix-A. 

We have studied the effect of temperature on the cor- 
relation functions XO) Xi and X2 and on the kinetic en- 
ergy. Fig. 10 shows the variation of the on-site correlation 
fimction (xo) and the kinetic energy (—{K}) with tem- 
perature for different g values. Except for strong cou- 
pling the Xo increases while the kinetic energy decreases 
with increasing temperature. It may be noted that when 
Xo increases, the values of xi and X2 decrease and the 
polaron size becomes smaller. For strong coupling the 
efi'ect of temperature is small but opposite. The above 
behavior point to the fact that the polaron gets more and 
more localized with increasing temperature in the regime 
of weak and intermediate e-ph coupling, while for very 
strong coupling the polaron size gets larger with increas- 
ing temperature. We have also studied the properties of 
the ground state and different excited states individually 
to get a clear understanding of the observed temperature 
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variation. It is found that in the regime of weak and in- 
termediate couphng the value of xa is larger while the 
kinetic energy is smaller for the excited states compared 
to those of the ground state and this leads to increase (de- 
crease) of xo with increasing temperature. We 
have given a representative plot showing the variation of 
Xo with energy for the ground and the excited states for 
3' = 1 in Fig. 11. For strong coupling [g > 2.5) we find 
that the xx> is smaller for the excited states (in general) 
than that for the ground state and the polaron becomes 
a larger polaron with increasing temperature. 

For a fixed temperature xa increases and the kinetic 
energy decreases with increasing g (Fig. 10), but their 
variation with g is much slower at higher temperatures. 
These results are qualitatively similar to that noted 
for the two-site system. It may be mentioned that 
Hohenadler^^ studied the 1-d Holstein model by quan- 
tum Monte Carlo method and observed similar variation 
of the kinetic energy with increasing temperature. 

In Fig. 12 we have shown the variation of Cy as a 
function of g for the ordered and disordered systems for 
T = 0.1. The specific heat shows a peak at intermediate 
coupling, which is suppressed with increasing disorder 
strength. The behavior of Cy for the site potentials (- 
0.5,0,0,0) is similar to that observed for the 2-sitc case 
with Cd = 0.5. In Fig. 12 we have also shown the plots 
of Xo"^ to show that the specific heat peak occurs in the 
region of g where Xo"^ ^-Iso undergoes a sharp change. 
However, the peak position in Cy, as mentioned previ- 
ously, would be mainly decided by the tuning of the en- 
ergy separation of two lowest levels of the system with 
the thermal energy at low temperatures. 

The kinetic energy contains contributions from both 
the coherent and incoherent hopping processes^". The 
Drude weight (in units of ne^) represents the coherent 
part of the kinetic energy. The contribution from the in- 
coherent hopping processes may be found out from the 
total kinetic energy and the Drude weight by using the 
f-sum rule (see appendix B), 



{K) {D) 



(15) 



where, (K) and (D) are the thermal average of the kinetic 
energy and the Drude weight, respectively, and 



^ Er,. - E„., 



(16) 



is proportional to the contribution to the kinetic energy 
from the incoherent hopping processes. 

In Fig. 13 we have shown the effect of disorder on the 
Drude weight of the four-site Holstein model for t = 0.5. 
The kinetic energy and the Drude weight are plotted 
against g for the ordered and disordered systems. At 
.9=0, the Drude wight and the kinetic energy have same 
values for the ordered case indicating that in absence of 
e-ph interaction and disorder the entire part of the kinetic 



energy comes from coherent hopping. In the range of in- 
termediate to strong coupling the Drude weight shows 
an exponential suppression. The kinetic energy, on the 
other hand, shows an exponential suppression only in the 
range of intermediate coupling and a 1 /g"^ behavior in the 
strong coupling region, as predicted by the strong cou- 
pling perturbation theory. For the disordered case the 
Drude weight is smaller than the kinetic energy even at 
g=0, because of disorder-induced incoherent hopping. 

We have investigated the effect of temperature on the 
coherent and incoherent parts of the kinetic energy by 
evaluating the kinetic energy and the Drude weight as a 
function of temperature for different values of g. Some 
representative plots (for g= 0.1, 0.5, 1 and 1.5) are given 
in Fig. 14. At low temperature the kinetic energy as well 
as the Drude weight show negligible dependence on tem- 
perature, represented by the flat region of the curves. 
This flat region is larger for smaller g. After the flat re- 
gion the kinetic energy reduces exponentially with tem- 
perature. At high temperatures, where the Drude weight 
is very small, the kinetic energy may be fltted to a func- 
tion a/T + bT. For small value of g the value of b is very 
small and the kinetic energy varies approximately as 1 /T. 
In the temperature range 1.5<A;sT<2.5 the values of a 
and b are, respectively, 0.477 and 0.003 for 5 = 0.1, 0.450 
and 0.0054 for g = 0.5, 0.376 and 0.012 for g = 1.0 and 
0.281 and 0.020 for g = 1.5. The Drude weight shows 
almost exponential dependence on temperature both in 
the intermediate and high temperature range. A good fit 
to the Drude weight is obtained with a function ae~^^ , 
but with different values of the parameters a and b in the 
two (intermediate and high temperature) ranges. 

The incoherent part of the kinetic energy can be di- 
rectly determined from the difference between the neg- 
ative kinetic energy and the Drude weight. In Fig. 14 
we have also shown the variation of the incoherent part 
as a function of temperature for different g values. Ex- 
cept for very strong coupling (not shown in the figure) 
the incoherent part increases with temperature, reaches 
a peak and then decreases at a slow rate with increasing 
temperature. A cross-over temperature (Tcross) niay be 
defined from the intersection of the curves for the coher- 
ent and the incoherent part of the kinetic energy such 
that for T < T^ross the coherent part is dominant and 
for T > Tcross the incoherent part is dominant. We have 
plotted the variation of Tcross with g in Fig. 15. The 
cross-over temperature decreases with g. The Tcross t^s. 
g curve shows a sudden change in the gradient at g=l 
for t = 0.5. The rate of fall of Tcross with g is higher in 
the region g > 1 than that for g < 1. 

To get the combined effect of the temperature and the 
disorder we have studied the Drude weight and incoher- 
ent part of the kinetic energy as a function of temper- 
ature for the site potentials (0,0,0,0) and (-1,0,0,0). In 
Fig. 16 we have given such a plot for intermediate cou- 
pling {g = 1) where the effect of disorder is large. It 
is seen that at low temperatures the effect of disorder 
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is very large on both the Drude weight and the inco- 
herent part of the kinetic energy. The Drude weight is 
rapidly suppressed while the incoherent part is enhanced 
a lot by the disorder site potential. At high temperatures 
where the Drude weight becomes very small, the disor- 
der has little effect on the incoherent part of the kinetic 
energy. The effect of the disorder on the kinetic energy 
is even smaller because in this range of temperature the 
Drude weight slightly decreases while the incoherent part 
slightly increases with the introduction of the disorder. 

We would address now the possible size effect on the 
results. In the Appendix-A we have derived a strong cou- 
pling effective polaronic hamiltonian where the different 
coefficients (in eq. A-2) have no size dependence pro- 
vided the number of nearest neighbors (z) is same. For 
a four-site system, z ~2 which is same as that for an in- 
finite 1-d chain. The kinetic energy in the ground state, 
obtained from Eq. (A2), for 2; = 2 is given by 



Kg 



-2tp - 4-^ 



E 



^2(n+m) 



LU ^ n\ m\ in + m) 

n,m=0;n+m>l 

(") 



n=l 



where = te~^ . The above clearly shows that for 

strong coupling the kinetic energy in the ground state 
should not have any size dependence if the number of 
nearest neighbors is same, hence the kinetic energy for 
= 4 will be same as that for an infinite lattice, where 
N is the number of lattice sites. This is completely con- 
sistent with the results obtained in a recent numerical 
study by Hohenberg et al.^^. They obtained the same 
kinetic energy for A'' = 4,8,16,32 for large g while a 
small size effect is observed for small and intermediate 
coupling. With increasing temperature the size effect is 
even smaller as noted in Ref.l9. 



increasing e-ph coupling, temperature and strength of 
the disorder. For small values of g and T this contri- 
bution (Drude weight) shows a weak suppression with g 
or T and then a rapid suppression with increasing g (or 
T). The contribution from the incoherent hopping, on 
the other hand, increases with increasing g ox T initially, 
reaches a maximum and then decreases. This contribu- 
tion decreases as l/g^ with g in the strong coupling. In 
the regimes of strong coupling or high temperature the 
coherent contribution to the kinetic energy is very small 
compared to the incoherent contribution and the latter 
becomes almost independent of the disorder strength. 
We have also identified a cross-over temperature below 
which coherent part of the kinetic energy is dominant 
and above which the incoherent part is dominant. This 
cross over temperature decreases with the increase of the 
electron-phonon coupling. 

The variation of the polaronic specific heat with re- 
spect to the e-ph coupling shows a peak in the inter- 
mediate e-ph coupling regime at low temperatures. A 
suppression of the specific heat due to disorder has been 
observed when e-ph coupling is in the weak or interme- 
diate regime. 

For the Holstein model the interaction is very short 
ranged and the size effect generally does not play a sig- 
nificant role in shaping different properties^^'^^'^^. As 
pointed out previously the size effect on the kinetic en- 
ergy is negligible for strong coupling, while the effect is 
finite but small for weak and intermediate coupling. At 
higher temperatures the size effect is even smaller. Re- 
garding the role of the disorder on the Drude weight of 
the polaronic system, the size effect has not been reported 
to our knowledge. We find that for strong coupling or at 
higher temperatures the effect of disorder on the polaron 
dynamics is negligible, hence the size effect will also be 
negligible as similar to that for the ordered system. In 
other region of the coupling the qualitative behavior of 
the results presented here for the four-site system are 
expected to be the same for larger systems. 



IV. Conclusions 

We have investigated the effect of disorder and tem- 
perature on the properties of a polaron for a two- and 
four-site Holstein model. It has been observed that both 
the disorder and temperature reduce the polaron size and 
suppress the kinetic energy in the weak and intermediate 
e-ph coupling regime. For strong e-ph coupling polarons 
are practically immobile and the kinetic energy arises 
mainly from the to and fro motion of the (bare) electron 
between nearest-neighbor sites keeping the centre of the 
polaron fixed. In this regime the effects of disorder and 
temperature on the correlation functions and the kinetic 
energy are very small. 

The polaronic kinetic energy has contributions from 
the coherent and the incoherent hopping processes. The 
contribution from the coherent hopping decreases with 



APPENDIX- A 

Applying the standard Lang Firsov (LF) transforma- 
tion to the Holstein model in Eq. (1) one obtains 

H = e^ife-^ 

+ uj^blbi (A-1) 

i 

where R = gni{bi — ^|), Cp = g^w is the polaron bind- 
2 

ing energy, tp = is the polaronic hopping strength 
and j is a nearest neighbor of i. 

The second term of Eq.(A-l) represents the kinetic 
energy operator which involves zero-phonon as well as 
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multiphonon processes (in the LF phonon basis) asso- 
ciated with the hopping of a polaron between nearest- 
neighbor sites. In the strong couphng Umit, where tp is 
very small, following a second order strong coupling per- 
turbation theory an effective expression for the kinetic 
terms (within the phonon ground state) may be obtained 



as 



23 



oo 



n,m— 0;n+m> 
^2m 



^ n! to! (n + to)^; 



_ 2^ V V 

^■^toItow-^ * 



(A-2) 



i,k 



where z is the number of nearest neighbors. The first 
term of Hkin represents the coherent hopping of polarons 
without emission or absorption of phonons, the second 
term originates from virtual hopping of a polaron from 
site i to a nearest neighbor site j and back. This hop- 
ping is associated with virtual emission of multi phonons 
at the sites i and j followed by absorption of all the 
emitted phonons while hopping back. In this process the 
(bare) electron within the polaron undergoes forward and 
backward motion between nearest-neighbor sites keeping 
the polaron immobile. Third term represents an effec- 
tive hopping of a polaron from site i to a second near- 
est neighbor site k via a common nearest neighbor site 
j where multi phonons are created and then absorbed. 
For strong coupling the coefScients of the first and third 
terms in Eq. (A-2) become negligible and only the sec- 
ond term (~ t^/cp) contributes to the kinetic energy. In 
presence of disorder (n -|- m)ui in the denominator of the 
second term modifies to {n + to)^; ± e^, where is the 
difference between the site potentials at sites i and j. 
The major contributions of this term for strong coupling 
comes from multi phonon processes with high values of 
{n + m)u> such that (n -I- m)uj >> td- As a consequence 
of the above reasons the coefficient of the second term in 
(A2) as well as the kinetic energy in the strong coupling 
region becomes very weakly dependent on the disorder 
potential. 

It may be mentioned that the Eq.(A2) is valid for any 
dimension and size of the system. The value of z would 
be different for different dimensions. 



APPENDIX- B 

The Hamiltonian of the system in the broken time re- 
versal symmetry may be written as 

H{(^) = - ^(te^-^cjci+i + h.c) + V (B-1) 



where V is composed of electronic potential and interac- 
tion of electrons with lattice. Let Hq be the Hamiltonian 
of the system when <^ = 0, and the corresponding eigen 
energies and eigen states be and |V'n)i respectively. If 
En{(f)) be the eigen energy of the Hamiltonian H{(j)), the 
Drude weight for the n-th eigen state of the Hamiltonian 
Ho is" 

^" - U=o 



For small (j), H{(j)) can be expressed as 
H{<p) = Ho- i(j)tJ2i4<^i+i - 4+i<^i) 

i 

+ ^t^{4ci+i + cl^iCi) + 



(B-3) 



Considering the terms involving as a perturbation a 
second order perturbation calculation gives 

Er^i^) = El- z0i(^O|(^clc,+i - 4iQ)|^0) 

i 

K€l(E.4cm-4+iCi)l€')P 



E°, - EO 

n' n 



+ ^(^^1 E(4ci+i + 4+iC,)|V°) + .... (B-4) 

i 

Therefore 



Dn = It" 



miEi4ci+l-cl,C,)\rn'}\' 

K'-K 



(B-5) 



where, 



Kn = I E(c1q+i + 4+,Ci)\i;l) (B-6) 



is the expectation value of the kinetic energy for the n-th 
eigen states of the Hamiltonian Hq. Taking the ther- 
mal average over all the eigenstates at a temperature T 
establish the sum rule 

2 2 



En — En' 



(B-8) 



where the thermal average of an observable A is 

^ ' En^M-En/kBT) ■ 

In the above expression is the expectation value of the 
observable A for the n-th eigen state of the Hamiltonian 
Ho. 
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FIG. 1. Convergence of the correlation function (niW2)/(ni) {u2 = 62 + 62) ^o'' two-site system, (a) Solid line: t= 2.1, 
ed= £2 — €1=1.0 with np=40; (b) dashed line : t=1.0, td= 0.5 with np=40; (c) dash-dotted line : t=1.0, ed=1.0 with np=40. 
Solid squares are for np=15, solid triangles are for np=20, solid circles are for np=25. All energy parameters are expressed in 
units of 0^=1. 
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FIG. 2. Convergence of the kinetic energy —{K)/t for the two-site system. Solid line: f=1.0, £^=1.0 with np=40, solid 
squares are for np=15, solid triangles are for np=20, solid circles are for np=25. All energy parameters are expressed in units 
of ui=l. 
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FIG. 3. Variation of {riiUj) /2g with g+ for t = 2.1 and ed=1.0 for the two-site Holstein model. Solid line:(niwi)/2g', dashed 
line:(n2W2)/2p and dash-dotted line:(niM2)/25=(n2Wi)/2g. 
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FIG. 4. Variation of the kinetic energy —{K)/t (solid lines) and the correlation function Xd"* (dashed lines) with g+ for 
t=2.1 and different for the two-site Holstein model. Solid lines from top to bottom and dashed lines from bottom to top are 
for £d.=0, 0.5, 1.0 and 2.0, respectively. 
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FIG. 5. Variation of Xd with gj^ for the ordered case (e^ = 0) for different temperatures. Solid line: ksT =0, long dashed 
line: ksT =0.5, short dashed line: fcsT =1.0, dotted line: fcsT =1.5, dash-dotted line ksT =2.0. t = 2.1 in the energy scale 

of 01=1. 
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FIG. 6. Variation of the kinetic energy —{K)/t with g+ for the ordered case for different temperatures. Solid line: fc^T =0, 
short dashed line: ksT =1.0, dotted line: ksT =1.5, dash-dotted line ksT =2.0. t = 2.1 in the energy scale of uj=l. 
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FIG. 7. Variation of the specific heat (Cij/fcs) with g+ for the ordered two-site Holstein model for different temperatures. 
Solid line: fcsT =0.1, long dashed line: ksT =0.2, short dashed line: ksT =0.3, dotted line: ksT =0.4, dash-dotted line 
feflT =0.5. t = 2.1 in the energy scale of uj=l. 
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FIG. 8. Variation of the specific heat (C'v/kB) with g+ at a temperature ^^^=0.1 for ordered and disordered two-site 
Holstein model for different hopping parameters. Solid line: t =1.0, £^=0; long dashed line: t =1.0,£d=0.5; short dashed line: 
t =2.1,ed=0; dotted line: t =2.1,6^=0.5; dash-dotted line t =2.1,6^=1.0. 
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FIG. 9. Variation of the kinetic energy —{K)/2t and the correlation function Xo"^ with g in the ground state for ordered and 
disordered four-site Holstein model. Solid linesiordered case, dashed lines: disordered with site potentials (1,0,0,0), dash-dotted 
curves: site potentials (-1,0,0,0). Solid triangles and diamonds represent xi and X2, respectively for the ordered case. Value of 
the hopping parameter used t = 1.0 in the energy scale of uj=l. 
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FIG. 10. Variation of xo and the kinetic energy —{K)/2t with temperature for the ordered four-site Holstein model for 
different values of g. Value of the hopping parameter used t = 1.0 in the energy scale of u!=l. 
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FIG. 11. Values of the on-site correlation function (xo) for the ground and the excited states for g = I and t = 1 for the 
ordered four-site Holstein model shown up to S = 2.4, where E is the energy of the excited state measured from the ground 
state energy. There are sixteen states within the energy range < E < 2.4. Four states have double degeneracy (shown by 
solid circles). 




FIG. 12. Variation of the specific heat (Cu/fcs) and the correlation functions Xo"^ with gr at a temperature fcBT=0.1 for 
t=1.2 for different set of site potentials for four-site Holstein model. Solid linesiordered case, dashed line: disorderd case with 
site potentials (-0.5,0,0,0) and dash-dotted line: disordered case with site potentials (-1,0,0,0). 
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FIG. 13. Variation of the kinetic energy —{K) /2t and the Drude weight with g in the ground state for ordered and disordered 
cases for t=0.5. Ordered case: Kinetic energy (sohd Une), Drude weight (dashed Une). Disordered case with site potentials 
(-1,0,0,0): Kinetic energy (dotted line), Drude weight (dash-dotted line). 




FIG. 14. Variation of the kinetic energy —{K)/2t (sohd lines) and the Drude weight {D/2t) (dashed lines) with ksT for the 
ordered lattice for different values of g. The dot-dashed lines represent the incoherent part of the kinetic energy. Value of the 
hopping parameter used t=0.5. 
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FIG. 15. Variation of the cross-over temperature (Tcross) with g for t=0.5 for the ordered case. 
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FIG. 16. Variation of the Drude weight and the incoherent part of the kinetic energy (in a scale of 2t) with temperature 
for ordered and disordered cases for t—0.5. Ordered case: Drudc weight (soUd hne), incoherent part of the K.E. (dashed hne). 
Disordered case with site potentials (-1,0,0,0): Drude weight (dotted line), incoherent part of the K.E. (dash-dotted line). 
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